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'We investigate the stability and phase transition of localized modes in Bose-Einstein Condensates 
(BECs) in an optical lattice with the discrete nonlinear Schrodinger model by considering both two- 
and three-body interactions. 'We find that there are three types of localized modes, bright discrete 
breather (DB), discrete kink (DK), and multi-breather (MUB). Moreover, both two- and three- 
body on-site repulsive interactions can stabilize DB, while on-site attractive three-body interactions 
destabilize it. There is a critical value for the three-body interaction with which both DK and MUB 
become the most stable ones. We give analytically the energy thresholds for the destabilization of 
localized states and find that they are unstable (stable) when the total energy of the system is higher 
(lower) than the thresholds. The stability and dynamics characters of DB and MUB are general for 
extended lattice systems. Our result is useful for the blocking, filtering, and transfer of the norm in 
nonlinear lattices for BECs with both two- and three-body interactions. 

PACS numbers: 03.75.Lm, 63.20.Pw, 67.85.De, 03.75.Hh 

Keywords: Bose-Einstein condensates, discrete localized modes, three-body interaction, stability 


I. INTRODUCTION 

Localized excitation in Bose-Einstein condensates (BEGs) has become one of the most interesting topics in nonlinear 
lattice systems since the discrete breathers (DBs) and the intrinsic localized modes are discovered [J-Q. The most 
well-known fascinating feature of the localized mode is that it can propagate without changing its shape as a result 
of the balance between nonlinearity and dispersion DB arises intrinsically from the interplay between the 

nonlinearity and the discreteness of the system. DB has been observed in vario us sy stems, such as micromechanical 
cantilever arrays Isl, antiferromagnet systems @,113: Josephson-junction arrays |lll fT3|. nonlinear waveguide arrays 
[H [g, BECs g g, Tonks gas g], superfluid fermi gases g|, and some dissipative system s g . g. The 
static, dynamical, and other properties of DB haye been studied theoretically in the last decade [2ll - [g . It has been 
demonstrated that DBs are attractors in dissipatiye systems [26l - ig . or act as virtual bottlenecks which slow down 
the relaxation processes in generic nonlinear lattices It is shown that the stability of the discrete localized 

modes plays a crucial role in blocking, filtering, and transfer of the norm through a localized mode. By far, there are 
many interesting works focused on this stability by considering two-body interactions 

Recently, the three-body interactions could be obseryed or realized in experiment and theory [H, [13. In 2014, 
Petroy [40|| proposed a method to control the two- and three-body interactions in ultracold Bose gas in any dimension. 
The three-body interactions play an important role in many interesting physical phenomena |4ll - l44l| . and eyen lead 
to a yariety of unique properties that are absent in the system dominated by the two-body interactions which can 
be governed by a Feshbach resonance [13. For example, in 2010, Dasgupta [13 discovered that if the two-body 
interactions are attractive, the presence of three-body interactions makes the crossover process from Bardeen-Gooper- 
Schrieffer (BGS) to Bose-Einstein condensates (BEGs) a nonreversible one. In 2012, Singh et al. [d^l found that the 
coupling of the two- and three-body interactions can affect strongly the transition from Mott insulator to superfLuid 
for ultracold bosonic atoms in an optical lattice or a superlattice. Up to now, there are few works on localized 
excitations in nonlinear lattice systems by considering three-body interactions. Especially, there are no systematical 
analysis of the types, existence, and stability of the localized modes, such as DB, the discrete kink (DK) [isl - ig . and 
multi-breather (MUB). It is natural to ask how the two- and three-body interactions affect these properties of the 
localized modes in BECs. 

In this paper, we investigate the stability and phase transition of localized modes in BECs in an optical lattice with 
a discrete nonlinear Schrodinger model (DNLS) in the case by considering both two- and three-body interactions. 
'We find that there are three different types of localized modes, that is, DB, DK, and MUB, and give the critical 
conditions for these localized modes. Both the two- and three-body on-site repulsive interactions can stabilize DB, 
while the three-body on-site attractive interactions destabilize it. We calculate analytically the energy thresholds, the 
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Peierls-Nabarro (PN) energy barrier characterizing the stability of the localized excitation modes. If the total 

energy of the system is higher (lower) than the thresholds, the localized states are unstable (stable). Moreover, the 
stability and dynamics characters of DB and MUB are general for extended lattice systems. Our result is important 
for the transfer of BECs through the discrete localized modes, and is useful for controlling the transmission of matter 
waves in interferometry and quantum-information processes when there are both two- and three-body interactions in 
the system. 


II. LOCALIZED STATES AND PEIERLS-NABARRO BARRIER 

A. The model 


Besides the on-site two-body interactions, let us investigate the effect on transfer of BECs through discrete localized 
mode from the on-site three-body interactions of ultracold Bose gas in an optical lattice. Under the mean-field theory, 
the Hamiltonian of the system of BECs in an optical lattice can be written as : 


^ JJ ^ TJ 


n—1 


n=l 


jM-1 

n—1 


( 1 ) 


Here n (= 1, • • • , M) is the index of the site, ipn is a complex variable and |'0n(t)P = Nn{t) is the mean number of 
bosons at site n (i.e., the norm Nn{t)). The first two terms represent the mean-field two- and three-body interaction 
energy, respectively, and the third term describes the hopping between nearest-neighboring sites. Ui = Anh’^agYeff /m 
represents the effective on-site inter-atomic two-body interaction, where 14// is the effective mode volume of each 
site, m is the atomic mass, and ag is the s-wave atomic scattering length. Here, we focus on the repulsive two-body 
interaction, i.e., Ui > 0. U 2 represents the effective on-site inter-atomic three-body interactions, including both the 
repulsive and the attractive three-body interactions which are represented by C /2 > 0 and U 2 < 0, respectively. J is 
the tunneling amplitude. Within the canonical equation one can obtain the dimensionless DNLS [SSl - l^ 

i^^=Xl ^i’n + X 2 IV'nl ^ -f , (2) 


where 


M 

E = 1- (3) 

n—1 

Here, Ai = Ui/J, X 2 = U 2 /J and t = Jt are the normalized dimensionless two-body interaction, three-body interac¬ 
tion and time, respectively. Assume that ijjnit) = A„(/)exp (z0„(t)), the Hamiltonian H becomes 

M . . M-l 

n^l^ n^l 

Usually, we use the Peierls-Nabarro (PN) energy landscape [5l|, to reflect the fact that discreteness breaks the 
continuous translational invariance of a continuum model. It is related to the PN potential whose amplitude can be 
seen as the minimum barrier which should be overcome to translate an object by one site. As in Ref. [s^, the PN 
energy landscape is defined as follows: for a given configuration of amplitudes A„, with respect to the phase difference 
50ij = Oi — 6 j, the PN energy landscape is obtained by extremizing H 

HpN = ^PN = max(-7J), (5) 

duij ouij 

where and Hp^^ are the lower and upper parts of the PN landscape, respectively. 

In order to give an insight into the dynamical behavior of BECs in an optical lattice with both two- and three-body 
interactions, we mainly consider the nonlinear trimer model, i.e., the DNLS with M = 3 lattice sites. In this case, 
the Hamiltonian H of the system is 

Hm=3 = {A\ -\- A 2 + Ag) -I- — (Ag -I- A® -I- A|) 

- [A 1 A 2 cos { 61 - d 2 ) + A 2 A 3 cos { 02 - ^ 3 )] , 


( 6 ) 
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where 60i2,S023 € [0,7r]. When 50i2 = 5023 = 0, one can get the upper PN energy landscape 

ifpN = {Af + ^2 + Af) —^ (A® + ^2 + ^ 3 ) 

+ (^ 1 +^ 3 ) ^ 2 - (7) 

When 5012 = 5023 = tt, the lower PN energy landscape can be obtained as 

ifpN = — — {A\ + A 2 + Ag) — (A® + ^ 1 ® + ^ 3 ) 

— {Ai+A3)A2. ( 8 ) 

The lower and the upper parts of the PN landscape bound the phase space of the trimer [s^. Because the localized 
mode whose properties we are studying corresponds to the minimum on we should focus on the lower PN 

landscape, i.e., 50i2 = 5023 = 


B. Phase transition of localized states and Peierls-Nabarro barrier 


To investigate the type of localized modes and its stability when the norm transfers through a localized mode, we 
use the nonlinear trimer model (M = 3 for Eq.([2|)) 

idtipi = Ai IV-il + A 2 '‘V'l - ^V'2, 

idti’2 = Ai |'02| '^i’2 + A 2 |V'2| '‘V'2 - ^ + i’ 3 ), (9) 

idtii3 = Ai \ip3\'^i’3 + A 2 IV'sl "‘V's - ^V'2- 


Here, the normalization reads 


^ = En=i IV'nl ^ = 1. By setting 50 


2 A 1 H 1 (1 - 2 A\ - Al) + 2 A 2 A 1 (1 - 


-^l-Al-Al 

2 A 1 H 3 (1 - 224^ - 


j Al (Hi + A 3 ) 
Al) + ‘2X2A3 (1 — 


-^l-Al-Al + 


A 3 {Al + 2 I 3 ) 

^^-Al- All 


12 = 5023 = from Eq. dH]) one can get 

Al){l - 2A! - Aj) 

= 0 , 

Ai){l - 2Al2 - Al) 

= 0 . 


( 10 ) 


Let us define 


and give inequality 


dAl 


= Hi, 


d Aid A 3 


= H2, 


dAl 


= H3, 


Hi - H 1 H 3 > 0. 


( 11 ) 


( 12 ) 


By solving Eq. OT one can get the different kinds of solutions, including the stationary states and the saddle points. 
By substituting the saddle points obtained from Eq. m into Eq. m, one can get the critical conditions which are 
the boundaries between phases I-III in Fig. 1. Actually, the boundaries can also be gained numerically from Eq. m- 
The results show that phase I occurs when Ai and A 2 satisfy the relation 


A 2 > 1.68708 - I.O 88 OIA 1 + 0.00522A?. 


(13) 


For this case the norm can be pinned in any one of the three sites, shown in Fig. 2(a)-(c). There exists DB. That is, 
the norm is pinned in the middle site. 

Phase III appears when Ai > 3.5 and Ai and A 2 satisfy the relation 


A; < A 2 < A„, 


(14) 
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FIG. 1: (Color online) The norm distribution and the structure of the localized states with two-body interactions (Ai) and 
three-body interactions (A 2 ). The solid lines are the boundaries dividing phases I-III. The dotted line presents the critical value 
A 2 with Ai. When A 2 = A 2 the DK is the most stable one. The symbols (a)-(h) correspond to the sets of parameters (Ai,A 2 ) 
used in Fig. 2. Phase IV is gained by Eq. 113 which corresponds to a four-site model. 


where 

A; = 7.39452 - 4.33978Ai -b 0.37649A? - 0.0I78A?, 

A„ = -12.223 -b 6.6531Ai - 1.722AJ -b 0.1715A^ 

In this case, the norms are localized in the adjacent two sites, shown in Fig. 2(e)-(g). It is called as DK. 

Phase II is the part other than phases I and III, shown in Fig. 1. For this case, the norms of the three sites are 
nearly equivalent and there is no saddle point on the contour plots of the lower PN energy landscape Hpj^, shown in 
Fig. 2(d) and (h). In this case, no localized mode exists. 

To investigate the stability of the localized states, we pay our attention to the PN barrier. As shown in Fig. 2, the 
projection of Hp]\j onto the Ai — A 3 plane exhibits one, two, or three minima. Each minimum refers to a stationary 
state. If there exist saddle points, the localized mode appears and these stationary states correspond to DB or DK. It is 
clear that the existence of DB, DK, and saddle points depends strongly on Ai and A 2 . We define the energy of DB (DK) 
asEoB {Edk), the energy of the saddle points as Etftr, and the energy difference AE = Ethr—EjjB(DK)- Actually, AA 
is the PN barrier. When the total energy of the trimer Et = —El > Ethr or Et — Ebb(dk) > Ethr — Ebb{dk) = AA, 
the DB (DK) should be unstable. On the contrary, when Et — Ep,B(DK) < the DB (DK) should be stable. 
We should note that Ethr merely marks the energy of the trimer at the saddle point and is identified with the 
destabilization threshold of the DB (DK). AE is the energy difference between the stationary state and the saddle 
points. Hence, AE represents the minimum energy barrier required to translate BECs by one lattice site. The higher 
AE is, the more stable the DB (DK) is. Therefore, the quantity AE can provide a deep insight into the stability 
property of DB (DK) for different values Ai and A 2 . 

The structure of MUB in the extended lattices is similar to that of a DK. MUB is a four-site solution of the DNLS, 
with high atomic density concentrated mainly at two middle sites and two low-density sites. Up to now, there are few 
studies on the properties of MUB, especially on its stability. Here we also investigate the effect on transfer of BECs 
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FIG. 2 : (Color online) Contour plots of the lower PN energy landscape for different A2 with fixed Ai = 6. (a)-(h) 

correspond to the points a-f marked in Fig. 1 . In (a)-(c), the three minima are separated by the two saddle points and DB 
exists. In (d) and (h), no saddle point exists, which means there is no any kind of localized modes. In (e)-(g), two minima are 
separated by one saddle point and DK exists. The color codes present the energy of Hpj^. 


through the MUB from the on-site two- and three-body interactions. It can be investigated by using four-site model, 
i.e., M = 4. Eq. ([9]) becomes 


idttpi = Ai + A 2 IV'il - ^"02, 

idtil2 = Al + A 2 |V'2| "^V'2 - ^ + ^ 3 ) , 

= Al |^/>3| V 3 + A 2 IV'sl ^ (^^>2 + ^a) , 

idtiIjA = Al I'tpil + A 2 \tpA\ "^V'4 - ^V'3- 


(16) 


Here, the normalization reads N = X]n=i IV'nl ^ = 1- Similarly, by setting 6612 = 6623 = S834 = tt, one can get 


\i[A^-A\)+X2{A\-A2) + -(^ 
Ai(a12 — ^3) + A 2 ( 4 l 2 — ^3) 2 ( 


1 / H2 ^1 + ^3 \ _ p, 
1 f Al + 743 A 2 + A 4 


A 2 


A?, 


= 0 , 


(17) 


Al( 4 l 3 — 414) -I-A2(H| — H4) 4 - - 


4 44\ , ( A 2 + Aj ^3 I _ Q 


A 3 


A, 


By solving Eq. (HZl, one can get the MUB and the saddle points. 

Now, let us discuss the stability and dynamics of DB and DK on phases I and III in which the system has saddle 
points (i.e., the system has the Peierls-Nabarro barriers), excluding phase II in which there is only one stationary 
state without saddle points. As shown in previous works [3, if there is no saddle point, the system should be in a 
random (generic) state in the presence of boundary or other local dissipation, and cannot form any kind of localized 
modes. That is, for the system with parameters of phase II, the localized mode does not occur. Although phase 
IV cannot be shown in Eig. 2 as there are not enough dimensions, saddle points can still exist and they are also 
investigated here. 
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FIG. 3: (Color online) Dynamics of the trimer when its total energy is increased with fixed two-body interactions Ai = 6 and 
three-body interactions A 2 = —1.5, where Ethri = —1.91469. A projection of the orbit onto the Ai — A 3 plane is over-plotted 
(black curve) in (a)-(d). (a) and (e) Et = —1.9176 < Ethri, areas in phase space are disconnected and the amplitudes An{t) 
indicate that the maximum amplitude remains at site 2, i.e., the DB is stable; (b) and (f) Et = —1.91179 > Ethri, areas in 
phase space are connected. A slight instability of the DB centered at site 2 is observed, the breather migrates to site 1, and 
then tangles in site 1; (c) and (g) for a high total energy Et = —1.8995, the norm is transmitted to site 1, tangles, and then 
comes back; (d) and (h) for an enough high total energy Et = —1.5504, the orbit explores large parts of the phase space and 
visits all three sites. In all cases 5e = tt. 


III. STABILITY OF LOCALIZED STATES 

The contours of the lower PN energy landscape Hpj^ for different A 2 with fixed Ai = 6 are shown in Fig. 2(a)-(h). 
Fig. 2(a)-(c) corresponds to phase I in Fig. 1 and there exist three stationary points and two saddle points. That 
is, the norm can be localized in three different ways (see Fig. 1) and DB exists. Fig. 2(e)-(g) corresponds to phase 
III in Fig. 1 and there exist two stationary points and one saddle point [this case cannot exist in the system without 
three-body interactions, as shown in Fig. 1 with A 2 = 0]. This case corresponds to DK. Fig. 2(d) and (h) corresponds 
to phase II in Fig. I and there exists only one stationary point but no saddle point. 


A. The stability and dynamics of DB 


To study the transfer of BEC through the DB in an optical lattice with three-body interactions, we should pay 
attention to the energy threshold Ethr and energy Edb again. In this case, for Ai —>■ 00 and A 2 /A 1 —)> 0, one can get 
the saddle point from Eq. (HOD as 


1 


^1 ~ 


A? 


A2 




V 2 XI 2^/2Xf ’ 
3 1 1 


(18) 


2y/2Xl V2Xl y/2Xi 2 V 2 XI ■ 


By substituting the saddle point into Eq. ([5D with Af + A 2 + = 1, one can get the energy threshold 


Etitri = + ^ 

4 2 4Ai 4Af 4Af 


16 At 


3 


1 

12 


1 


1 


3 5 

4 A? SAl'^SA?/^ 


16A? 16Af 2Af 


Al + 
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32Af 4Af 


Ai 


(19) 
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FIG. 4: (Color online) (a) Energy difference AE and (b) perturbations SAi for destabilizing the DB as a function of three-body 
interactions A 2 . Different lines indicate three different values of Ai, respectively. The left and right of the dotted line present 
attractive (i.e., A 2 < 0) and repulsive (i.e., A 2 > 0) three-body interactions, respectively. S9i2 = n/4 in (b). 


Of course, Edb can also be obtained by substituting the bright breather given by Eq. m into Eq. ([ 8 ]). 

Next, we consider the fixed point corresponding to the bright breather which is gained from Eq. m- An initial 
condition for the bright breather reads , A^^). We add perturbations to site 1: -ip (t = 0) = 

((Af’'®-|-( 5 ^Je*'^'’,A 2 ,Af’®), where A 2 = (1 —IV'iP — Compared to the bright breather, we add an amplitude 

5ai to site 1 and rotate the phase 61 by 5e. Dynamics on the PN landscape for increasing total energy of the trimer with 
fixed two-body interactions Ai = 6 and three-body interactions A 2 = —1.5 is shown in Fig. 3. Here Ethri = —1.91469, 
and we fix = tt and increase 5ai in Fig. 3(a)-(d). If the perturbation is small, Et < Ethri, the areas in phase space 
are disconnected. Furthermore, the DB is stable and practically no transfer of norm takes place on short time scales 
[see Fig. 3(a) and (e)]. If the perturbation is large, Et > Ethri, the areas in phase space are connected [see Fig. 3(b) 
and (c)]. Instability of the DB centered at site 2 can be observed, that is, the breather migrates to site 1 and norm 
is transferred to site 3 [see Fig. 3(b), (c), (f), and (g)]. If the perturbation is large enough, the orbit goes out of 
the regular island into the chaotic sea [see Fig. 2(d)] and large amplitudes A„(t) are found at all the three sites, as 
depicted in Fig. 3(h). 

Further, we use the parameters 6 Ai and AE to investigate the effect of the two- and three-body interactions on the 
stability of DB, shown in Fig. 4. It is clear that both AE and 6 Ai increase with Ai whether three-body interactions 
are repulsive or attractive (i.e., A 2 > or A 2 < 0), which means that a larger perturbation is required to destabilize the 
DB when Ai increases and on-site two-body interactions can stabilize the DB. Interestingly, both AE and 5Ai increase 
with repulsive three-body interactions (A 2 > 0 ) but decrease with attractive three-body interactions (A 2 < 0 ), which 
means that a relatively large perturbation is needed to destabilize the DB when repulsive three-body interactions 
increase. When attractive three-body interactions increase, a relatively small perturbation is needed to destabilize 
the DB. That is, repulsive on-site three-body interactions can stabilize the DB, while attractive on-site three-body 
interactions destabilize the DB. For large enough attractive on-site three-body interactions, AE = 6 Ai = 0 and the 
DB is completely unstable. 


B. The stability and dynamics of DK 


One can investigate the stability and dynamics of DK with the same way used above. When A 2 —>■ —00 and 
A 1 /A 2 —>■ 0, one can get the saddle point from Eq. (fTUl) as 


1 Vs 54V3 


3V3Ai ^ 513V3 
SXj 256A^ ■ 


Ai — A3 


V 3 8A2 128 A 2 


( 20 ) 











FIG. 5: (Color online) Dynamics of the trimer when its total energy is increased with fixed two-body interactions Ai = 6 and 
three-body interactions A 2 = —7.5, where Ethr 2 = —1.44281. A projection of the orbit onto the Ai — A 3 plane is over-plotted 
(black curve) in (a)-(d). (a) A contour plot of the lower PN energy landscape Hpj^ is shown for the total energy of the trimer 
below the threshold {Et = —1.44304 < Ethr 2 )- Obviously, the areas in phase space are disconnected, (e) The amplitudes An{t) 
with time t indicate that the dominating amplitude remains at sites 2 and 3, i.e., the norm is still localized in sites 2 and 3 and 
the DK is stable, (b) and (f) for total energy Et = —1.44238 > Ethr 2 , areas in phase space are connected. A slight instability 
of the DK centered at site 1 is observed. The DK migrates to site 1 and then tangled at site 1 or 3, but the norm of site 
2 is nearly constant, (c) and (g) for a high total energy Et = —1.43809 > Ethr 2 , the DK migrates to site 1 and then come 
back more easily. Long-range and long-lived Josephson oscillations between sites 1 and 3 with negligible variation of norm in 
site 2 are observed, i.e., the DK is unstable and the dominating amplitude does not remain at sites 2 and 3. (d) and (h) for 
an enough high total energy Et = —0.871835, the orbit is out of the regular island into the chaotic sea, and the dominating 
amplitude can be found in any of the three sites, which means the DK is unstable. In all cases Se = 0. 




FIG. 6: (a) (Color online) Energy difference AE and (b) perturbations 5Ai for destabilizing the DK as a function of three-body 
interactions A 2 . Different lines indicate three different values of Ai, respectively. A 2 is a critical value of A 2 . When A 2 = A 2 , 
the DK is the most stable one for the fixed Ai. 59 i 2 = 0 in (b). 


By substituting the saddle point into Eq. ([5]) with A\ 4 






29 


one can get the energy threshold 
981 151119 


512A^ 


16384A^ 


1 3 45 9369 3013119^ 

6 16A| 32A3 1024A4 32768A5 ) ^ 

135 243^^2 ( 27 189 


16A| 
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Next, we consider the fixed point corresponding to the DK which is also gained from Eq. OT. The initial 
condition for the DK reads = 0) = , A?^). One can add perturbations to site 1: V'(0) = 

{{Ai^ + , ^ 2 , ^ 3 *^), where A 2 = {1 — Dynamics on the PN landscape for increasing total 

energy of the trimer with fixed two-body interactions Ai = 6 and three-body interactions A 2 = —7.5 is shown in Fig. 
5. Here Ethr 2 = —1.44281. If the perturbation is small, Et < Ethr, the areas in phase space are disconnected. The 
DK is stable and the norm is localized nearly in the adjacent two sites 2 and 3 on short time scales, shown in Fig. 
5(a) and (e). If Et is just larger than the areas in phase space are connected and instability of the DK can be 
observed; that is, the DK migrates from site 3 to site 1 and then tangles at 1 or 3, but the norm of site 2 [presented 
by A 2 {t)] is nearly constant, shown in Fig. 5(b) and (f). If the perturbation becomes larger, one can see that there 
are long-range and long-lived Josephson oscillations between sites 1 and 3 with negligible variation of norm in site 2, 
shown in Fig. 5(c) and (g). If the perturbation is large enough, the orbit explores large parts of space and visits all 
three sites, shown in Fig. 5(d), and large amplitudes An{t) can be found at all three sites, shown in Fig. 5(h). 

One can also use SAi and AE to investigate the stability of DK, shown in Fig. 6. For a fixed Ai, one can see that 
both AE and SAi first increase and then decrease with A 2 . That is, there exists a critical three-body interaction A 2 
with which the DK is the most stable one. When A 2 A 2 , the DK becomes more unstable when A 2 decreases or 
increases. Furthermore, we numerically get the critical value when Ai > 3.5 

a; = 0.27034-I.I 3 I 86 A 1 -b0.00516A?. (22) 

It is presented by the dotted line in phase III in Fig. I. 

Obviously, the stability of the DK depends strongly on both two- and three-body interactions, while the DK is 
more unstable than the ordinary DB. 


C. The stability and dynamics of MUB 


For MUB, the lower PN energy landscape in Eq. o becomes 


iJpN— ^(^1+^2+^3 + ^4) —^ i^i+ ^2 + M+^T) 

— {^AiA2 -\- ^ 2^3 + ^ 3 ^ 4 ) • 


(23) 


Ai 
A 3 
Ai 

By substituting the saddle point into Eq. 


1 get the saddle point from Eq. ( 

fT71) as 

1 

73 
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373 Ai , 

51373 


8 A 2 

128Ai 

-U 

8 Ai 

256Ai ’ 

1 

73 

5473 

373 Ai , 

51373 

73 + 

8 A 2 

128Ai 

SAi 

256Ai ’ 

27^3 





4Ai 





. dSSl) 

with 

A! + Ai 

+ Al + Al = 

= 1 , one 

^2 

2 

1 Ai 

403 -f 12Ai 

27 

3 + 

8 A 2 6 

64Ai 


3 (-8541 - 

f 3072Ai 

-b 640Ai) 


2048Ai 

9 (l03827-b61338Ai-b8032Af-3584A?) 


(24) 


8192At 


(25) 


Suppose that the initial condition for the MUB is = 0) = (A^^^, A^^^, A^^^, One can add 

perturbations to site 1: ^ {0) = {{A’^^^ + , A 2 , A^’^^, where H 2 = (I — — |^/’ 3 p — |^/> 4 p)^/^. 

Dynamics of the MUB for increasing total energy of the four-site system with Ai = 6 and A 2 = —7 is shown in Fig. 7. 
Here Ethrs = —1.50027. If there is no perturbation, i.e., SAi = S6 = 0, the MUB is stable absolutely and the norm 
is localized in the middle adjacent two sites 2 and 3 with A 2 {t) = A 3 {t) and Ai{t) = Ai^t), shown in Fig. 7(a). If 
the perturbation is small, Et < Ethrs, the MUB is stable and the norm is localized nearly in the adjacent two middle 
sites 2 and 3, shown in Fig. 7(b). In contrast, if Et is just larger than Ethrs, the MUB becomes unstable after t ^ 78 
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FIG. 7: (Color online) Dynamics of MUB for increasing the total energy of the four-site system with Ai = 6 and A 2 = —7, 
where Ethrs = —1.50027. (a) 5a = Se = 0, the MUB is absolutely stable, and A\{t) = A 4 ,{t) and A 2 {t) = Ai{t). (b) Sg — -k jZA, 
Et = —1.50212 < Ethr 3 , the multi-breather is stable, and it is still located at sites 2 and 3. (c) Sg = 7r/2.8, Et = —1.47925 is 
just larger than EthrS, the multi-breather is unstable, and the norm from sites 2 and 3 migrate to sites 3 and 4 after t « 78 time 
steps, (d) Sg = 7r/2, Et = —1.42046 > EthrS, the stability of the multi-breather is destroyed, and the dominating amplitude 
can be found in any of the four sites. In all cases 5 ai = 0. Initial conditions are given by Eq. (I17II . 



FIG. 8: (Color online) (a) Energy difference AE and (b) perturbations SAi for destabilizing the MUB as a function of three- 
body interactions A 2 . Different lines indicate three different values of Ai, respectively. The A 2 is a critical value of A 2 . When 
A 2 = A 2 , the MUB is the most stable for the fixed Ai. 59 i2 = 7r/3.4 in (b). 


time steps, and the MUB migrates from sites 2 and 3 to sites 3 and 4 or 1 and 2, shown Fig. 7(c). Similarly, if the 
perturbation is large enough, the stability of the MUB is destroyed, shown in Fig. 7(d). 

Also, one can use SAi and AE to investigate the stability of MUB, shown in Fig. 8. For a fixed Ai, both AE and 
SAi first increase with A 2 before decreasing. That is, there still exists a critical three-body interaction A 2 at which 
the MUB is the most stable one. When A 2 A 2 , the MUB becomes unstable when A 2 decreases or increases. 


IV. EXTENDED LATTICES 


Now, let us generalize our study to the case with extended lattices, i.e., M > 4. 
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FIG. 9: (Color online) The stability and migration of a DB in extended lattices (M = 101 sites). The color code represents 
\ipn{t)\. The insert in (a) represents the structure of DB. Initial conditions are given by Eq. (I26II with (a) 5a = Se = 0, 
Et = —2.61029 < Ethri, (b) Sa = 0.184, Se = -k, Et = —1.91179 is just larger than Ethri, (c) Sa = 0.199, 5g = tt , 
Et = —1.87136 > Ethri- In all cases Ai = 6 and A 2 = —1.5. 


A. DB in extended lattices 


Here we use the same initial conditions as the case with three sites, i.e., Ai = 6, A 2 = —1.5, to study an extended 
lattice system. We assume M = 101 and the DB locates at the site n = 51. The initial condition reads 


V^5o(0) = (Af^+<5AJe*^^ 

■*/'5l(0) = A 2 , 

1/152(0) = 

ipn(0) = 0, else. 


(26) 


where A 2 = (1 — ll/'soP — and Af^ is obtained exactly from Eq. (nnj. The wave function is normalized 

to ll/'riP = 1- The stability and migration of a DB in extended lattices (M = 101 sites) is shown in Fig. 9. 

The solution for the DB in the trimer, i.e., the initial condition shown in Eg. (1261) . is inserted in the extended lattice. 
Here Ethri = —1.91469 when Ai = 6 and A 2 = —1.5 (corresponding to Fig. 3). Obviously, when no perturbation is 
added to the site 50, i.e., 5ai = = 0, the breather is stable, shown in Fig. 9(a) and (d). When the perturbation is 
added to the site 50 and the total energy of the local trimer Et is just larger than Ethri, the breather is stable and 
no migration takes place, different from the case with three sites, shown in Fig. 9(b) and (e). The reason is that the 
energy can flow into additional degrees of freedom in extended lattices. When the perturbation is large and the total 
energy of the local trimer Et > Ethri, the breather is destabilized and migrates from site 51 to site 50 after t « 2 
time steps, similar to the case with the three-site model, shown in Fig. 9(c) and (f). 
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FIG. 10: (Color online) The stability and migration of a MUB in extended lattices (M = 101 sites). The color code represents 
\tpn{t)\. The insert in (a) represents the structure of MUB. Initial conditions are given by Eq. H17|l with (a) Sa = Se = 0, 
Et = —1.55596 < Ethrs, (b) Se = tt/2.8, Et = —1.47925 is just larger than Ethri, (c) 5g = 7r/1.79, Et = —1.39563 > Ethri- In 
all cases Ai = 6 and A 2 = —7. 


B. MUB in extended lattices 

For MUB, we use the same initial condition as the case with the four-site model to study its stability and dynamics 
in the extended lattice system (M=101 sites), i.e., Ai = 6 and A 2 = —7. Here, the MUB is located at the sites n = 51 
and n = 52. The initial condition reads 


■*/'5l(0) = A 2 , 

^(0)=^™, (27) 

^53(0) 

V’n(O) = 0, else, 

where A 2 = (1 — IV’soP — |l/' 52 p — and and A^^^ are obtained exactly from Eq. (flTl) . The wave 

function is normalized to IV’nP = 1- The solution for the MUB in the four-site model is inserted in the extended 

lattice. Here Ethrs = —1.50027 when Ai = 6 and A 2 = —7. Obviously, when no perturbation is added to site 50, 
i.e., 5ax = de = 0, the breather is stable, shown in Fig. 10(a) and (d). As |V'5o(OI = \'^ 53 {t)\ and \ip 5 i{t)\ = |V' 52 (t)|, 
the red and the blue lines are overlapped, and the black and the green lines are overlapped in Fig. 10(d). When the 
perturbation is added to site 50 and the total energy Et is just larger than Ethrs, the energy can flow into additional 
degrees of freedom in extended lattices, and the MUB is still stable and no migration takes place, shown in Fig. 10(b) 
and (e). When the perturbation is larger and the total energy Et > Ethrs, the MUB is destabilized and migrates from 
site 51 to site 53 and locate in sites 52 and 53 after t k, lA time steps, shown in Fig. 10(c) and (f). 

From the discussion above, one can see that the stability and the dynamics characters of both DB and MUB are 
general for extended lattices. 
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V. DISCUSSION AND SUMMARY 

Although the three-body interactions could be observed or realized in experiment and theory [H, , there are 

no systematical analysis of their type, existence, and stability of the localized modes. Previous works indicated that 
DB can exist and its stability plays a crucial role in blocking, filtering, and transfer in BECs [29l - l3Ti [s^ . The exact 
condition to destabilize DB is still not clear, especially in the case by considering three-body interactions. Recent 
works [HHi showed that MUB can exist for BECs in optical lattices by only considering two-body interactions in 
the case with large interactions and appropriate phases. 

In our work, we systematically investigated the existence of the different localized modes for different parameters 
(Ai, A 2 ), i.e., DB, DK, and MUB, and discussed the effect of two- and three-body interactions on their stabilities. Both 
on-site repulsive two- and three-body interactions can stabilize DB, while on-site attractive three-body interactions 
destabilize DB. DK and MUB are the most stable ones when the three-body interactions are equal to a critical value 
for the fixed two-body interactions. It may provide effective guidance to gain different kinds of localized modes in 
optical lattices BUI' and help us to study its properties in experiment, especially for MUB. 

Besides, in our work, three analyzed thresholds to destabilize the three localized modes are given explicitly. If the 
total energy of the system is higher than the energy thresholds, the localized state is unstable. On the contrary, 
if the total energy of the system is lower than the energy thresholds, the localized state is stable. It may lead 
to some interesting applications for blocking and filtering atom beams when there are both two- and three-body 
interactions in the system. Furthermore, it is useful for controlling the transmission of matter waves in interferometry 
and quantum-information processes 

In summary, we have investigated the stability and phase transition of localized modes in BECs in an optical lattice 
with the discrete nonlinear Schrddinger model by considering two- and three-body interactions. It has been shown 
that there are three different types of localized modes. The first one is bright DB which can be stabilized by both 
on-site repulsive two- and three-body interactions. However, on-site attractive three-body interactions destabilize 
DB. The second one is DK which is the most stable one when the three-body interaction is equal to a critical value 
for fixed two-body interactions. The third one is MUB. It becomes the most stable one when three-body interaction 
is in the critical value. Moreover, the stability and dynamics characters of DB and MUB are general for extended 
lattice systems. Our results provide a deep insight into the dynamics of blocking, filtering, and transfer of the norm 
in nonlinear lattices for BECs by considering both two- and three-body interactions. 
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